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We present fits to charmless hadronic B decay data from tlie BaBar, Belle and Cleo experiments 
using two theoretical models: (i) the QCD factorization model of Beneke et al. and (ii) QCD fac- 
torization complemented with the so-called charming penguins of Ciuchini et al. When we include 
the data from pseudoscalar-vector decays the results favor the incorporation of the charming pen- 
gum terms. We also present fit results for the unitarity triangle parameters and the CP-violating 
asymmetries. 
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I. INTRODUCTION 

A wealth of experimental data on hadronic charmless B decays has become available from the BaBar and Belle 
experiments. These studies of the numerous B decay channels are designed to test the Cabbibo-Kobayashi-Maskawa 
(CKM) explanation of CP violation in the standard model as represented by the unitarity triangle condition VudV*f^ + 
^cdV*ij + VtdV^l = on the CKM mixing parameters. Such tests have been facilitated by the recent significant 
progress in the theoretical understanding of hadronic decay amplitudes based upon QCD factorization which allows 
the amplitudes to be expressed in terms of numerous soft QCD parameters, such as meson decay constants and 
transition form factors, and a set of calculable coefficients. In this paper we present an analysis of recent data, 
based upon QCD factorization. We also investigate the potential contribution to the decay amplitudes of b quark 
annihilation and so-called charming penguins. The data that we attempt to fit includes decays to pseudoscalar (tt and 
K) and vector [p,Lj,K* and (f>) mesons. This is an extension of an earlier study Q| that was based upon simplified 
formulae derived from the heavy quark limit of QCD factorization. 

The starting point for the calculation of all B meson decay amplitudes is the effective Hamiltonian 

Heff(/^) = ^ ^ Ap[Ci(Ai)gf + C2(M)Q?+ C,{^l)Q,]+ other terms (1) 

V ^ p=u,c j:=3,...,10 

where Xp = V*^Vpb is a product of CKM matrix elements, q = d,s and the local AB — 1 four-quark operators are 

Ql = {qaPa)v~A{Ppbi3)v~A, 
Ql = {qaPl3)v-A{Pl3ba)v-A, 
Q3,5 = iQf3bl3)v-AY{qWa)vTA, 

q' 

Qi,Q = {qpba)v-A y^(gagg)vTA; 

q' 

Q7,9 = {qi3bp)v-AY{q'aq'a)v±A: 

q' 

(98,10 = {qi3ba)v~AY{qW0)v±A, (2) 

where q' G {u,d,s,c}, a and /3 are color indices, eq> — 2/3(— 1/3) for u{d)-type quarks and we use the notation, for 
example, 

Q5 - {<iphp)v^A J2^qWJv+A = [ga7p(l " I5)ba] [q'pl^ {I + Ts)?^ • (3) 
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Qi.2 are the current-current tree operators, Q3....,6 are QCD penguin operators and Q7....,io are electroweak penguin 
operators. The 'other terms' indicated in |^ include the electromagnetic and chromomagnetic dipole transition 
operators. In the Standard Model the contributions from the electroweak penguin operators and magnetic dipole 
operators are generally small. The exception is the electroweak penguin operator coefficient Cg which is larger than 
the QCD penguin coefhcients C3 and C5. 

It has been shown by Beneke et al. 2] that, in the heavy quark limit, the hadronic matrix elements of the four-quark 
operators in the amplitudes for non-leptonic B decays into two light mesons A/1.2 have the form 



and can be calculated from first principles, including non-factorizable strong interaction corrections. QCD factorization 
extends naive factorization by separating matrix elements into short distance contributions at scale 0{\/m\,) that are 
perturbatively calculable and long distance contributions 0(1/Aqcd) that are parameterized. 

B meson decay can also be initiated by b quark annihilation with its partner. Although the annihilation contributions 
to the decay amplitude are formally of ©(AqcD/^^fc) and power suppressed, they violate QCD factorization because 
of end point divergences. However these weak annihilation contributions can be included into the decay amplitudes 
by treating the end point divergences as phenomenological parameters. Analyses based upon QCD factorization with 
inclusion of weak annihilation have been undertaken for B PP HQ and B PV Although general agreement 
with experiment was found, some branching ratios for PV decays were only marginally consistent and all predictions 
were plagued by the large uncertainties associated with poorly determined parameters within the theory. A recent 
global analysis of PP and PV decays |5j found that QCD factorization plus weak annihilation could fit many decay 
channels but yielded results too low for the B nK* decays. More recently, Aleksan et al. Q have undertaken 
a global analysis of PV decays and have concluded that QCD factorization cannot fit the experimental data when 
decay channels involving K* mesons are included, and suggest that this failure is due to some larger than expected 
nonperturbative contribution. Motivated by the concept of so-called charming penguins, that is non-perturbative 
©(Aqcd/™?)) corrections from enhanced c-loop penguins, first introduced by Ciuchini et al. they introduce 

additional long-distance contributions to the decay amplitudes and include the two complex parameters from these 
additional amplitudes in their global fit. They obtain a slightly better fit but their best-fit parameters are at the limits 
of the allowed domain. A recent study 0, limited to i? — > p7r decays, has used QCD factorization to place bounds 
on the Fi^'^ form factor. A detailed study of QCD factorization applied to PP and PV decays of B mesons has 
just been completed by Beneke and Neubert Predicted branching ratios and CP asymmetries for a large number 
of PP and PV channels are given for default values of input parameters and detailed estimates of the theoretical 
uncertainties in these predictions determined for various scenarios of input parameters. Beneke and Neubert find 
that there is a scenario for which there is general global agreement between the results of QCD factorization and 

measurement except for B tt^K and the group of B ^ ttK* decays. 

In this paper we undertake a global analysis of 18 PP, PV and VV channels using two theoretical models, (i) QCD 
factorization with inclusion of weak annihilation and (ii) this model supplemented with charming penguins. Our study 
is similar in spirit to that of Aleksan et al. for PV decays but we extend the global fit to include PP and some 
VV decays. This paper is organized as follows. In Sec. II we review the structure of the decay amplitude within QCD 
factorization and discuss the various parameters that occur in this amplitude. Inclusion of weak annihilation and 
charming penguins is discussed in Sec. Ill and Sec. IV respectively. The method and results of our best fit for our 
two models to current experimental data is presented in Sec V, and Sec. VI contains our discussion and conclusions. 



In QCD factorization, the amplitude for B decay into two light hadrons (mesons) Af 1 2 has the form, neglecting 
weak annihilation processes. 



{M,M2mB) = (Afi|Ji^|S)(A'f2|J2''|0) l + Y 



rna'^ + 0{AQCD/mt) 



(4) 



n 



II. DECAY AMPLITUDE IN QCD FACTORIZATION 



(A//iM2|Heff|S) 



7f 



{ ^ A„ar [T^ (Ml ,M2)+ T, {M2 , Afi )] 



i=l,2 




(5) 



p—u,c z— 3,...,6,9 



where 



TM,M2) = (A//i|jl7'^(l-75)5|S)(Af2|g7M(l-75)u|0), 
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TABLE I: Leading order Wilson coeSicients d in the NDR scheme calculated at the scales ^ and fj,h ~ \/AhfJ., where = 0.5 
GeV, for n — mt and fi = m(,/2. The input parameters are J^qqI^^'' — 0.225 GeV, mt{mt) = 167.0 GeV, mb{mb) = 4.2 GeV, 
Mw = 80.42 GeV, a = 1/129, and sin^ 9w = 0.23. 



Scale (GeV) 


4.20 


2.10 


1.45 


1.02 


Ci 


1.1174 


1.1848 


1.2392 


1.3120 


C2 


-0.2678 


-0.3873 


-0.4755 


-0.5862 


Ca 


0.0121 


0.0185 


0.0235 


0.0299 


C4 


-0.0274 


-0.0383 


-0.0459 


-0.0551 


Cs 


0.0080 


0.0105 


0.0120 


0.0136 


Ce 


-0.0341 


-0.0526 


-0.0677 


-0.0883 


Cr/a 


-0.0140 


-0.0282 


-0.0314 


-0.0303 


Cg/a 


0.0288 


0.0432 


0.0555 


0.0734 


Cg/a 


-1.2913 


-1.3785 


-1.4408 


-1.5190 


Cio/a 


0.2888 


0.4213 


0.5194 


0.6424 



TABLE II: Next to leadini 


g order Wilson coefficients d 


in the NDR scheme calculated at the scales 


fj, and = V^hMj where 


Ah = 0.5 GeV, for /x = rub 


and fi = mb/2. The input parameters are Aq^^'' 


= 0.225 GeV, mt(mt) = 


167.0 GeV, mb{mb) = 4.2 


GeV, Mw = 80.42 GeV, a 


= 1/129, and sin^ 6*^ = 0.23. 








Scale (GeV) 


4.20 


2.10 


1.45 


1.02 


Ci 


1.0813 


1.1374 


1.1820 


1.2405 


Ca 


-0.1903 


-0.2948 


-0.3700 


-0.4619 


Ca 


0.0137 


0.0212 


0.0274 


0.0358 


C4 


-0.0357 


-0.0506 


-0.0618 


-0.0762 


Cs 


0.0087 


0.0102 


0.0105 


0.0096 


Ce 


-0.0419 


-0.0653 


-0.0854 


-0.1146 


Cr/a 


-0.0026 


-0.0139 


-0.0147 


-0.0100 


Cs/a 


0.0618 


0.0986 


0.1317 


0.1813 


Cg/a 


-1.2423 


-1.3181 


-1.3526 


-1.4149 


Cio/a 


0.2283 


0.3388 


0.4015 


0.4991 



T2(Mi,M2) = (Afi|g7'^(l-75)6|B)(M2|M7^(l-75)u|0), 

T3(Mi,M2) = (Mi|g7^(l-75)6|B)(M2|g'7M(l-75)9'|0), 

T4(Mi,M2) = (Mi|gV(l-75)^>|S)(M2|g7M(l-75)9'|0), 

T,iM,,M2) - (Mi|g7^(l-75)&|S>(M2k-'7M(l + 75V|0), 

Te{M^,M2) = -2(Afi|g'(l-75)6|B)(M2|g(l+75)<z'|0), 

n{M^,M2) = (Mi|g7''(l-75)&|S)(M2|e,,g'7^(l-75)9'|0). (6) 

The two-quark matrix elements are the well determined electroweak decay constants f-mfK^fp, e.t.c. and the B 
transition form factors F^^, Fk^ Ap, e.t.c. In principle transition form factors are independently measurable in B 
semi-lcptonic decays but to date they are only loosely constrained by measurements and model estimations. 

The coefficients a.i have to be calculated from the Wilson coefficients '2\ . We have calculated the Wilson coefficients 
at several scales of 0{mh) using the next-to-leading-order (NLO) renormalization group equations 

C,(m) = U,,{fi, Mw)Cj{Mw). (7) 

We follow the Beneke et al. [2| prescription of (i) dropping terms of 0{a1), O(a^) and O(asa) in (|7|), (ii) neglecting 
the effect of the electromagnetic penguins C7....sq{Mw) on the evolution of the QCD penguin coefficients Ci,...,6 and, 
(iii) in Cj{Mw), splitting the 0{a) electroweak penguin terms into those enhanced by large mt or 1/sin^ O^r^ which 
are treated as leading order (LO), and treating the remainder, together with the Oias) terms, as NLO. Our calculated 
values are shown in tables and ^ and are very similar to those of Beneke et al. and Du et al. Q . 

To lowest order in the strong coupling constant as, the coefficients are the same as in naive factorization, that is 

aLO = C, + aV7Ve (8) 

where i' = i — (—1)* and A'';; = 3 is the number of quark colors. The higher order corrections include lowest order 
gluon exchange between the quarks in the basic tree amplitudes which are calculated and folded into the light cone 
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TABLE III: 


Factorization ai,i{M2 


) coefficients evaluated at the scale fi using the expressions of |3,la 


|. The light cone corrections 


a^'^ are set 


equal to zero, making 


all coefficients except ae,/ universal. 








A on 


z.iu 


0-1,1 




1 or 70 1 A 0000^' 

i.Uo lZ + U.UzUUl 


i.U (91 + U.UoDyi 


0-2,1 




o nof;o n OQQf;^' 
U.UUOU — U.Uoobz 


n nQ'7'7 nil Qo^' 


03,1 




n nnc^c 1 n nnoi,>' 
U.UUOO + U.UUzil 


U.UUoo + U.UUoU 


Oi,I 




— V.UolZ — U.UiOii 


— U.Uoof5 — U.UzUo? 


Oa,i 




— U.UOD9 — U.UUuOi 


— U.U4io — U.UU/yz 






—0 nr)7n — n nn9fi7 


—0 ni 06 — n oo50j 


ohiP) 




-0.0433 - 0.0152i 


-0.0586 - 0.0188i 


ol,i{P) 




-0.0465 - 0.0056i 


-0.0630 - 0.0056i 


oliiV) 




-0.0075 - 0.0007i 


-0.0094 - 0.0013i 


ohiV) 




0.0009 - 0.0115i 


0.0019 + 0.0152i 


09,1 




-0.0094 - 0.0002i 


-0.0097 - 0.0003i 



TABLE IV: Factorization ai,j{Al2) coefficients evaluated at the scale fj, = 2.10 GcV using the expressions of 0,0 and Lii and 
Yang ilOl values for the light cone corrections a^^| . 





TT 


K 


p, CJ 


K' 


a 1,1 


1.0809 + 0.0369i 


1.0762 + 0.0432i 


1.0798 + 0.0369i 


1.0753 + 0.0439i 


"■2,1 


-0.0432 - O.liaOi 


-0.0290 - 0.1322i 


-0.0400 - 0.1130i 


-0.0261 - 0.1344i 


as,! 


0.0085 + 0.0037i 


0.0080 + 0.0043i 


0.0084 + 0.0037i 


0.0079 + 0.0043i 


"■4,1 


-0.0306 - 0.0206i 


-0.0313 - 0.0210i 


-0.0325 - 0.0206i 


-0.0321 - 0.0210i 


"■2,1 


-0.0343 - 0.0072i 


-0.0364 - 0.0061i 


-0.0386 - 0.0076i 


-0.0383 - 0.0061i 


0-5,1 


-0.0108 - 0.0050i 


-0.0102 - 0.0059i 


-0.0107 - 0.0050i 


-0.0101 - 0.0060i 


"■6, 1 


-0.0586 - 0.0188i 


-0.0586 - 0.0188i 


-0.0097 - 0.0012i 


-0.0102 - 0.0013i 


"■6,1 


-0.0630 - 0.0058i 


-0.0630 - 0.0058i 


0.0022 + 0.0147i 


0.0032 + 0.0152i 


09,/ 


-0.0097 - 0.0003i 


-0.0097 - 0.0004i 


-0.0097 - 0.00031 


-0.0097 - 0.0004i 



distribution functions ^M{x,fJ.) of the participating mesons (see, for example, 0). This results in the Ui coefficients 
having the form 

a,(AfiM2) = a,j{M2) + auiiMiN^) (9) 

where Mi is the recoil meson containing the spectator (anti) quark and M2 is the emitted meson. The complex 
quantities aij describe the formation of M2, including nonfactorizable corrections from hard gluon exchange or light 
quark loops in penguins. They do not involve the hard gluon exchanges with the spectator quark, these are described 
by the (possibly) complex quantities j/. 

In the correction terms the leading-twist light cone distribution functions for both pseudoscalar and vector mesons 
are expanded in the first few terms of a Gegenbauer expansion 

$m(2:,m) =6a:(l-a:)[l + ^afC3/2(2x-l)]. (10) 

n 

The asymptotic limit <i>7\/(a;, /i) = 6a;(l — x) is valid for the mass scale /i ^ cx). The parameters a^^ are anticipated to 
be small but they are not well established. To economize in the number of fitting parameters in the initial fits to data 
presented here they are taken to be zero. With this simplification all light mesons included in our analysis have the 
same spatial wavefunction and all coefRcients ai,i except ag,/ are the same for all decays. Formulae for the evaluation 
of the tti,! can be found in Although Beneke and Neubert obtain a different expression for ae,i{M2 = V) 

to that of Du et al. 4], this is not important here as aQ{PV) does not occur in the decay amplitudes for M2 = V since 
{y\{<il)s±p\^) = 0. The 0{as) corrections to a\'~' include contributions from one-loop vertex corrections Vm(^,q;*^) 
to all ai,i and from penguin corrections PM{fJ-i , Sp), involving the quark mass ratio Sp = {nip/mb)'^, to a4,i and 
oe,/. We neglect the small electroweak penguin corrections to Pm. Typical parton off-shellness in the loop diagrams 
contributing to Vm and Pm is 0{mb), suggesting that the scale /i (x rrib should be used in evaluating ai,i. The results 
of our calculations of a^./ coefRcients without light cone corrections are given in table IIIll Results including the light 
cone corrections af^j taken from 10] are given in table II V I The coefficients most affected are a2,i and 04./. 

The coefficients j/ are not universal even when light cone corrections are neglected. They contain not only 
the low energy parameters (decay constants and form factors) from the lowest order calculations, in common with 
naive factorization, but low energy contributions to the folding integrals involving another non-perturbative complex 
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parameter Xh which is only loosely constrained by model estimations. To discuss the form of the a^j/, we first note 
that, from the contributions of the coefficients to the decay amplitude for B M1M2 are of the form 



2 F i i 



(11) 



where g\ 2 and g\ are products of Clebsch-Gordan coefficients tablulated in, for example, P, US- Using the 
formulae of 0, we can write, with Nc — i, 



/Ma^Miaj,// — -z-^iCi>asl3i 



(12) 



where = +l(i — 1, 



,4, 9), £5 



-1, eg = 0, and 

fsfAhfMi 



3(1 + eia^ 



M2 



(1 



.M2 



..M2 



M2 



)(1 



Ml 



Ml 



The chiral factors r*^^ are zero for Mi a vector meson and are 



2mi 



„K 



2m\ 



mb{mu + rus) 



(13) 



(14) 



for the pseudoscalar mesons. It should be noted that these Uiji contributions to the decay amplitudes are independent 
of the B transition form factors. However, they do involve the poorly determined parameter fs/^B where Jb is the 
B leptonic decay constant and As ~ 0.6 GeV is related to the B light cone distribution function. The parameter 
X^^^ is the contribution of a logarithmic end-point divergence in the integration over the Mi light cone distribution 
function 



- 



dx 



1 



(15) 



These functions take no account of the internal quark transverse momenta which, if included, would make the integrals 
finite but not calculable within pcrturbative QCD. X^^ is parameterized as 



X'J^ = In 



niB \ 



A 



QCD 



(16) 



where pn is not expected to be larger than 3. We take ln(ms/AQCD) = 3.03. The energies involved in the calculation 
of Uiji imply that the appropriate scale is not that of the scale p, used in calculating the aij but ph = V ^hP where 
A/i = 0.5 GeV. We use this with p = mi,/2 in our fitting so that UsfB/i'mB^B) = 0.0209. We note that substantial 
light cone corrections a*^2 can significantly enhance the aiji coefficients. 



III. ANNIHILATION CONTRIBUTIONS 



Because of the heavy b quark mass it is expected that perturbative QCD calculations will give a reliable estimate of 
the annihilation contribution to the decay amplitude. In these calculations the basic perturbative quark amplitudes 
are again folded into the participating meson light cone distribution functions. Apart from the low energy regions of 
the folding integrals the only low energy parameters that appear in the lowest order calculations are the participating 
meson electroweak decay constants fB,f-K,fp, e.t.c. Again the low energy contributions to the integral introduce 
another non-perturbative complex parameter Xa, which is only loosely constrained by the model estimations. Detailed 
formulae are to be found in |2 0| . In this paper we follow the more extensive calculations of Du et al. ^ but express 
the annihilation contribution to the decay amplitude in the form 

{MiM2\niT\B) = ^BmiM, {A„[diCi + d2C2Wi - \t{di[C3^i + 

+ [diCA + deCeJAl + dz[C^ + N,C&]Ai} (17) 
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TABLE V; 


Annihilation coefficients di{Mi, M2) for B - 


-* M1M2. 


The VV channels refer to 


zero helicity states only. 




M1M2 


di 


d2 


di 






"6 


+ - 
TT TT 


-1 





— 1 


—2 


— 1 


—2 


n IT 


-1 





-1 

— i 


— Z 


1 

— 1 




TT TT 








U 


U 


U 


U 




1 





1 
X 


9 
z 


1 

— 1 


9 


p 7r+ 


1 





1 
X 


9 


1 
1 


9 


p TT 


1 





1 
1 


9 


fl 
U 


9 
















-^/2 

V 







u 


U 















U 


V ^ 


V2 













1 





1 


2 


-1 


2 



























I/V2 


1/V2 





-1/V2 





ujK° 








1/V2 





-1/V2 





n°K- 












-1/V2 





u 


— i 


-1 





-1 










1 


1 





1 





TT i\ 










U 




U 


P ^ 





I/V2 




U 




U 


P 





l/\/2 




n 


1 /,/9 


u 


■K+R- 








-1 





-1 





-K+K*- 








1 





1 





p+K~ 








1 





-1 





(jiK- 





1 


1 





1 





(f)K*- 





1 


1 





-1 













1 





1 













1 





-1 






where 

Cf 

BM1M2 = -TTa/s/Afi/Ma (18) 

c 

and Cf = {N^ — ^)/2Nc. The quantities A^{^^^{Ml, M2), where the superscript i{f) denotes gluon emission from initial 
(final) state quarks, result from folding the quark amplitudes into the meson distribution functions. If the asymptotic 
form of iCni) is used then 







18 (x.- 4+^) 








(r^^-r^^)(Xl-2XA + y), 


4{Pl.P2) 


= 67ra,(r^i +r^^)(2Xl -Xa), 


Ai(p,y) 




18 X,-4+!L 




A^(F,y) 


= vrasT 


^[2^2 + 6(X1-2Xa)], 


4{P.v) 


— dnas 


r^(2Xl-XA), 






= A\{P,V), 






= Ai{V,,V2) = 0. 





The coefficients di{Mi, M2) are Clebsch-Gordan type factors and are given in table^for the particle sign conventions 
Note that, in using (|17|l with table IVl there is no need to distinguish between PV and VP decays. 
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IV. CHARMING PENGUINS 



In our attempts to fit the data set within the scheme outlined above we found that tlic PP branching ratios could 
be accomodated with acceptable values of the CKM parameters and transition form factors. However, data on the 
ttK* channels is consistently too large to be accounted for. The problem is that the ttK branching ratios are only 
marginally larger than the ttK* ratios. In the QCD factorization scheme described above the penguin operators Q4 
and Qq contribute coherently and almost equally and dominate the ttK decay amplitudes whereas Qg is missing from 
the amplitude for irK* decay. This results in the predicted ratio of the ttK and ttK* branching fractions being too 
small. Perhaps, staying within the QCD factorization scheme, this failing can be removed by taking radically different 
light cone distribution functions for the K and K* mesons. We investigate the possibility of significant additional 
contributions from so-called charming penguins |l2l |. 

The largest term in the effective Hamiltonian that produces a strange quark comes from 

H = -J=V*sVcb[Cl{c0bi3)v-A{SaCa)v-A + C2iSaba)v ~Aic0Cp)v -a] ■ (20) 

The charming penguins originate from these terms when the c and c quarks annihilate. In the general description of 
two body B decays as given by Buras and Silvestrini charming penguins have the topologies CPi and DP2 of 
connected and disconnected penguins respectively. Ciuchini et al. 12] consider the contribution of these terms to tttt 
and 'kK decays in the 5L/(3) limit. In their notation, and including the small contribution from the u quark loop, 
this results in a contribution to the decay amplitudes which they express as 



A.. = -ml[VMP;+VubV:^lf^]g, (21) 



wherePi and pp^^^ are two complex numbers that are independent of the particular channel, tttt or nK. The only 
channel dependence is through the Clebsch-Gordan factor g which is the same as the Clebsch-Gordan factor in the 
04 contribution from QCD factorization. Ciuchini et al. suggest that all chirally suppressed terms should be dropped 
and replaced with this term. 

We take the charming penguin contribution to be from the penguin topology but to be in addition to the QCD 
factorization of the penguin. However, to retain the notation of 0, we express Pi and P^^^ as 

Pi = ^{UF^)\D\e^^, 

PGIM ^ 9^l^f^F^)\DGTMy'PG,M ^ (22) 

v2 



With the factor {frrFTr), taken here to be 0.042 GeV, the dimcnsionless parameters \D\ and must be less than 

unity as the charming penguins are expected to be small O^Aqcd/'tib) corrections. 

This simple model must be extended to include vector mesons. To this end we note that (i) in PV decays the 
vector meson must have zero helicity and that, for example, the s and q quarks forming the decay meson M can be 
expected to have zero spin projection in their direction of motion irrespective of whether they form a pseudoscalar 
or vector meson and (ii) that, when folded into the same light cone distribution functions, the amplitudes would be 
the same. This most simple model extends the SU{3) symmetry to SU{6). With this albeit simple extension, the 
charming penguin contribution to a particular amplitude is obtained from the factorization contribution by reference 

to the aij term. For example, the decay amplitude for B — s- ir^ p~ contains a term 



G 
72 

The charming penguin contribution is obtained by replacing this with 



l^f^F,[V^,V:,{alj - alj) + Vt^V;, alj]. (23) 



- < [V^bV:^ Pf + Pi] . (24) 

In the charming penguin amplitude, the s quark is produced from a left-handed field and can be expected to have 
predominantly negative helicity. The qq pair emanates from either right-handed or left-handed fields and, with zero 
helicity for the produced (sq) meson, we expect that the left-handed contribution will dominate to form an 04 type 
term. The right-handed term will be suppressed by a factor Aqcd /"ti b ■ We appreciate that the expectation of 
considerable suppression is false for the corresponding factorization term qqQq for which the suppression factor 
is of order unity. QqQq is a product of matrix elements of local operators containing {M\(sq)\0) which is zero for M 
being a vector meson. We suspect that the chiral enhancement of scalar meson production in factorization penguins 
is a property of factorization of the local product rather than a general feature of all charming penguin contributions. 
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TABLE VI: Measured branching ratios Br(exp), experimental error a, best fit theoretical branching ratios Br(BBNS) and 
Br(CP) for the Beneke et al. model (BBNS) and charming penguin model (CP) respectively, and the contribution to for 
various B decay channels. All branching ratios are in units of 10"''. For the VV channels the predictions are for longitudinal 
polarization states only as these decays are expected to be dominant. 



Decay 


Br(cxp) 


cr 


Br(BBNS) 


2 

X 


Br(CP) 


2 

X 


T 

TT TT 


4.8 


0.5 


5.0 


0.2 


4.9 


0.0 




5.6 


0.9 


4.0 


2.0 


5.9 


0.0 




25.4 


4.2 


24.4 


0.1 


23.3 


0.2 




16 


5 


7.3 


2.7 


11.1 


0.7 


P TT 


9.4 


2.0 


9.5 


0.0 


9.6 


0.0 




6.4 


1.3 


7.4 


0.5 


6.4 


0.0 




12.9 


1.2 


13.0 


0.0 


12.9 


0.0 


TT K 


18.2 


1.7 


20.9 


2.1 


19.5 


0.5 


IT K 


12.4 


2.6 


4.4 


7.7 


9.1 


0.8 


ujK- 


3.1 


1.0 


3.4 


0.1 


5.0 


2.2 


<j>K- 


8.8 


1.1 


8.4 


0.1 


8.4 


0.1 


(j>K'— 


5.4 


2.4 


8.9 


1.0 


9.0 


1.6 


TT+K- 


18.5 


1.0 


18.7 


0.1 


19.0 


0.2 


TT+K'- 


16 


6 


4.1 


3.8 


9.5 


1.0 




10.4 


1.4 


7.4 


3.9 


7.1 


4.6 




6.5 


1.7 


2.4 


5.2 


4.4 


1.2 




8.4 


1.6 


7.8 


0.2 


7.8 


0.2 


<fK'° 


7.2 


1.7 


8.3 


0.3 


8.3 


0.2 



V. FITTING METHOD 

We have attempted to fit the theoretical expressions for branching ratios with the available data as averaged by 
the Heavy Flavour Averaging Group ,13]. Measured branching ratios for 18 channels are shown in table IVll We take 
the measured branching ratios to be the mean of the B and B decays. For the two vector-vector channels (f)K* we 
multiply the measured branching ratios by the longitudinal polarization as measured by BaBar |l4l | . CP asymmetries 
are not included in the fit. The measurements are not always consistent between the experiments and the errors are 
large. We therefore prefer to compare the measured results with the predictions from the fit. To economize in the 
number of soft QCD parameters we have not included decay channels involving 77 and ?/ mesons. These amplitudes 
involve the mixing angle between the (uu + dd) and ss combinations. Also, in principle, there is mixing with cc 
which, though small, could make a significant contribution to decay modes through the enhanced quark decay modes 
6 cqc. 

For convenience we assign to each channel (Mi, M2) a number a. The statistical and systematic errors have been 
combined p!^ into a single error ctq.. The systematic errors are in general small and we ignore any correlations. We 
then form a function 

X^iPi) = '^[\Bra{Pi) - Bra{exp)\/aa]'^ + additional constraints. (25) 

a 

Bra{Pi) are the theoretical branching ratios expressed in terms of ten parameters Pi which we take to be the three 
Wolfenstein CKM parameters {A, p, 77} and seven soft QCD parameters {rj, r^, i^^, Fr-, ^p, ^a'* }• For the fit 

to the charming penguin model we introduce four more parameters: the modulus and phase of Pi and P^p^*^. In 
this case we fix A to the world average of 0.82 and keep rj'^ fixed at 1.0. The well established decay parameters 
{Itt, fx, fp, fuj, f(/), fK*} are held at their mean values and the Wolfenstein CKM parameter A is set to 0.2205. The 
results are not very sensitive to the divergence parameters and we held them fixed at pH = 2e'' ''* and Xa — 1.85e^'*^', 
values suggested by some preliminary investigations. Additional terms were included in the to take into account 
experimental and theoretical constraints from outside the data on B decay branching ratios. We search for a minimum 
of as a function of the Pi using the MINUIT |16|| program. 

Next to the experimental error on the measurement we have to consider the error from our assumptions on the QCD 
parameters that we do not fit for. Table IVTTl shows their central value and an estimate of the allowed variation. First 
we performed the fit using experimental errors only. With the best value we calculated a set of reference branching 
fractions. We then varied the parameters according to Table fVTll while keeping the value of the CKM parameters p 
and r] fixed. For each parameter this leads to a difference for each branching ratio. For every branching ratio we sum 
the difference in quadrature and consider this to be the model uncertainty. We add this uncertainty in quadrature to 
the experimental error and repeat the fit. 

The theoretical branching ratios and the contributions of the individual channels to based upon these best fit 
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TABLE VII: Value of the parameters used with their variation. 


Parameter 


Central Value 


Variation 




0.1307 


±0.00046 


fli 


0.1598 


±0.00184 


fp 


0.216 


±0.005 




0.194 


±0.004 


fK' 


0.216 


±0.010 


u 


0.233 


±0.004 


X 


0.2205 


±0.0010 




1.081 


±0.015 




0.5mi, 


mb 


\Ph\ 


2.0 


±1.0 


arg(pH) 


4.7 


1.6 




1.85 


±1.0 


arg(XA) 


2.86 


3.7 


' X 


1.0 


±0.2(CP only) 


Awoli 


0.82 


±0.05(CP only) 


„M 
Qn 





Table IV (Lii and Yang) 



TABLE VIII: Best fit values and one standard deviation errors for the Beneke et al. (BBNS) and charming penguin (CP) 
models. -FV./c and Ap^^^K* are the transition form factors F^^^{Q) for P and V mesons respectively, and rj'^ are the chiral 
enhancement factors which are nominally power suppressed but are Oil) in practice. 



Model 




Fk 


Ap 




A„, 


BBNS 
CP 


0.244 ± 0.038 
0.291 ± 0.022 


0.369 ± 0.031 
0.349 ± 0.077 


0.344 ± 0.098 
0.320 ± 0.065 


0.300 ± 0.094 
0.298 ± 0.072 


0.321 ± 0.136 
0.309 ± 0.117 






r'^ 


A 


P 


V 


BBNS 
CP 


1.09 ± 0.24 
1.0 


1.24 ± 0.17 
1.0 


0.813 ± 0.045 
0.82 


0.068 ± 0.071 
-0.044 ± 0.112 


0.383 ± 0.090 
0.397 ± 0.050 




\D\ 


Arg(D) 








CP 


0.068 ± 0.007 


1.32 ± 0.10 


0.32 ± 0.14 


1.00 ± 0.27 





values are given in table IVll The best fit parameter values are shown in table IVIIll together with our estimates of the 
errors. These errors are of course highly correlated. Plots of the error matrix ellipse for the Wolfenstein parameters 
p and fj are shown in figures [T] and □ The CKM angles are a = (78 ± 9)°, = (22 ± 2)°, 7 = (80 ± 7)° for the 
BBNS model and a = (63 ± 7)°, /3 = (21 ± 2)°, 7 = (96 ± 6)° for the CP model. These can be compared with the 
world averages of a = (96 ± 13)°, f3 = (23.3 ± 1.5)°, 7 = (61 ± 12)°. All errors are one standard deviation. For 
both models the results in table IVIIII for the values of the various form factors lie within the spread of theoretical 
estimates. The x^/dof is 34.5/17 for the BBNS fit and 14.5/13 for the charming penguins. 



VI. DISCUSSION AND CONCLUSIONS 



The first conclusion is that the factorization approach works quite well. Most of the branching ratios in table IVTI 
are predicted correctly by both models. The fitted parameters in table IVIIll look reasonable for both fits. The for 
the charming penguin fit is significantly better, due mainly to the poor fit of the BBNS model for the decay modes 
involving the K* meson. Also, the experimental value for cuK'^ is not easily accomodated within the BBNS model. 
Both models have a problem in fitting the tt'^K'^ mode. Figures and |21 show the position of the apex of the unitarity 
triangle for both fits. The results are consistent with each other for the angle 7 but give a larger value than that 
suggested by the CKM Fitting Group. The angle f3 agrees well for both fits. 

Our theoretical best-fit values for those CP asymmetries that have been measured are shown in table IIXI The 
direct CP-violating parameter Acp for the decay channel M1M2 is 

r{B ^ M1M2) +T{B ^ M1M2) ^ ' 

where B = bu 01 bd. The definitions of the other CP- violating parameters can be found in j^. Regarding these 
asymmetries it is too early to reach a conclusion. In many cases the different experiments disagree, in others the 
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TABLE IX: Measurements and theoretical best fit values for CP asymmetries. Only statistical errors are shown. 



^CP 


BaBar 


Belle 


BBNS 


CP 


K+-K- 


-0.102 ±0.05 


-0.07 ±0.06 


0.00 


-0.08 


K + 7V° 


-0.09 ± 0.09 


0.23 ±0.11 


0.06 


-0.02 


K°n+ 


-0.17 ±0.10 


0.07 ±0.09 


0.01 


0.11 


7r+7r° 


-0.03 ±0.18 


-0.14 ±0.24 


0.00 


0.00 


p+7r~ 


-0.22 ± 0.08 




-0.03 


-0.03 


p+K- 


0.28 ±0.17 




0.10 


-0.40 




-0.30 ±0.25 


-0.77 ±0.27 


0.01 


-0.23 


q 




—1 9^ ± n 41 


n 9n 


n 9^1 




0.36 ±0.15 




0.03 


-0.31 




0.19 ±0.24 




0.38 


0.65 


ACp^ 


0.28 ±0.19 




0.08 


0.31 




0.15 ±0.25 




0.02 


0.11 




-0.80 ± 0.38 


0.56 ±0.41 


-0.01 


-0.19 




-0.18 ±0.51 


-0.73 ± 0.64 


0.73 


0.62 


TABLE X: 


Predicted branching ratios (in units of 10 ^) and CP asymmetries for some channels not included 


in the fit. For 


the VV channels the predictions are for longitudinal polarization states only 


as these decays are expected to be dominant. The 


predictions 


are for the central values of the model fits with la errors estimated by sampling the parameter error matrix. 




Mean Branching 


; Ratio 


Acp 




Decay 


BBNS 


CP 


BBNS 


CP 




0.5 ±0.1 


1.1 ±0.2 


0.59 ±0.08 


-0.63 ±0.15 


P°P° 


0.5 ±0.1 


1.2 ±0.3 


0.69 ±0.16 


-0.66 ± 0.07 


n°p- 


8.4 ±0.7 


11.9 ± 1.4 


0.04 ±0.01 


0.06 ±0.01 


P°P- 


19.0 ±3.1 


17.8 ± 3.3 


0.00 ± 0.00 


0.00 ±0.00 


p°K+ 


2.0 ±0.6 


4.6 ±0.5 


0.00 ± 0.09 


-0.01 ± 0.64 


p°K*+ 


4.6 ±1.3 


5.2 ±0.7 


0.33 ± 0.06 


-0.35 ±0.06 



errors are so large that a meaningful discrimination is not possible. The theoretical asymmetries are very sensitive to 
the parameters and to the different models and, with improved statistics, could become the final test of factorization. 
In table we show our predictions for branching ratios and CP asymmetry Aqp for some channels not included 
in the fit ||22]. We include some VV channels that are currently under investigation. For these channels we have 
assumed that the decays are to longitudinally polarized vector mesons as it is expected 0, that decays to the other 
polarization states will be suppressed by at least a factor of (my/ms)^. 

In their most recent work Beneke and Neubert @| give formulae for the weak annihilation functions A^'3(P, V) that, 
in addition to the terms given in (|19|1 . include terms containing a parameter r^. This parameter has a similar 
origin to and, like r^, is suppressed by a power of Aqcd/'ti;, but, unlike r^, is not chirally enhanced. We only 
became aware of j8| after completion of this study. To check for the effects of these terms we have modified our 
program to include the contribution to annihilation and have reminimized x^- For the charming penguin model 
we found the best fit occurred with Xa — 1.09 exp(2.79i) and that there were very small changes in the best fit 
parameters of table IVIIII and the results of tables IVII IIXI and H For the BBNS model we found similar results but 
the overall fit was better in that the was reduced from 34.5 to 27.5. 

Finally, the results presented here are slightly different from preliminary numbers presented earlier |18| . for three 
reasons. The Heavy Flavour Averaging Group has updated some of the results and included a new branching fraction 
(p^K^). We have also corrected the vector- vector channels iC*^ for the effect of polarization. Finally we now include 
the systematic uncertainties in the fit. 



[1] W. N. Cottingham, L B. Whittingham, N. de Groot and F. Wilson, J. Phys. G:Nucl. Part. Phys. 28, 2843 (2002). 
[2] M. Beneke, G. Buchalla, M. Neubert and C. T. Sachrajda, Nuc. Phys. B 606, 245 (2001). 
[3] D. Du, H. Gong, Y. Sun, D. Yang and G. Zhu, Phys. Rev. D 65, 074001 (2002). 



11 
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7r"7r° = (1.9 ± 0.5) x 10"" and B" p"7r° = (11.0 ± 2.7) x 10"®. The former is on the high side of our predictions, 
particularly for the BBNS model, the latter is in good agreement with the prediction from the CP model but less so with 
that of the BBNS model 



